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Abstract. Here we model a star forming factory in which the continuous creation of stars results in a highly 
concentrated, massive (globular cluster-like) stellar system. We show that under very general conditions a large- 
scale gravitational instability in the ISM, which triggers the collapse of a massive cloud, leads with the aid of 
a spontaneous first generation of massive stars, to a standing, small-radius, cold and dense shell. Eventually, 
as more of the collapsing matter is processed and incorporated, the shell becomes gravitationally unstable and 
begins to fragment, allowing the formation of new stars, while keeping its location. This is due to a detailed 
balance established between the ram pressure from the collapsing cloud which, together with the gravitational 
force exerted on the shell by the forming cluster, acts against the mechanical energy deposited by the collection 
of new stars. We present a full analysis of feedback and show how the standing shell copes with the increasing 
mechanical energy generated by an increasing star-formation rate. The latter also leads to a rapidly growing 
number of ionizing photons, and we show that these manage to ionize only the inner skin of the standing star- 
forming shell. We analyze the mass spectrum of fragments that result from the continuous fragmentation of the 
standing shell and show that its shape is well approximated at the high mass end by a power law with slope -2.25, 
very close to the value that fits the universal IMF. Furthermore, it presents a maximum near to one solar mass 
and a rapid change towards a much flatter slope for smaller fragments. The self-contamination resultant from the 
continuous generation of stars is shown to lead to a large metal spread in massive (~ 10^ Mq) clusters, while 
clusters with a mass similar to 10^ Mq or smaller, simply reflect the initial metalicity of the collapsing cloud. 
This is in good agreement with the data available for globular clusters in the Galaxy. Other observables such as 
the expected IR luminosity and the equivalent width caused by the forming clusters are also calculated. 

Key words. Stars: formation; superstar clusters; supernovae: general; ISM: bubbles; globular clusters: general; 
Galaxies: starburst 



1. Introduction 

There is a mode of star formation that leads to massive 
(10^ M0 - a few 10^ Mq), densely concentrated collec- 
tions of stars. These have been named young massive clus- 
ters, super-star clusters and for the more massive ones 
the term starburst has also been used. The clusters are 
believed to have evolved from an interstellar cloud mass 
distribution N oc M^^ (Elmegreen & Efremov 1997) and 
thus young clusters present a similar power law (see also 
Zhang & Fall 1999) independent of the environment, while 
old (globular) clusters have a mass distribution that falls 
off towards low masses (Harris & Pudritz 1994), perhaps 
due to evaporation within a Hubble time (Elmegreen & 
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Efremov 1997). As pointed out by Ho (1997) young super- 
star clusters are overwhelmingly luminous concentrations 
of stars that present a typical half-light radius of about 
3 pc, and a mass that ranges from a few times 10"* Mq 
to 10^ Mq. The brightest ones have luminosities up to 
two orders of magnitude higher than R136 in 30 Doradus. 
Similar super-star cluster properties have been inferred 
from HST-STIS observations of AGN (Colina et al. 2002), 
and from radio continuum measurements of ultra-compact 
HH regions not visible in optical images, which indicates 
that they are the youngest, densest and most highly ob- 
scured star formation events ever found (Kobulnicky & 
Johnson 1999; Johnson et al. 2001). The massive concen- 
trations imply a high efficiency of star formation which 
even after long evolutionary times permits the tight con- 
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figuration that characterizes them, despite the impact of 
photo-ionization, winds and supernovae, beUeved to effi- 
ciently disperse the gas left over from star formation. Thus 
the self-gravity that results from the high efficiency of star 
formation is what keeps the sources bound together. This 
high efficiency also brings us to believe that the forma- 
tion of young clusters should be either a delayed or a very 
rapid event, to avoid negative feedback (Larsen & Richtler 
2000). The observational evidence now points to such mas- 
sive units of star formation (~ — 10^ M©) present 
at the excitation centers of blue compact and starburst 
galaxies such as M82 (de Grijs et al. 2001, O'Connell et 
al. 1995), and NGC 253 (Watson et al. 1996) as well as 
in galaxies of different types (see also Larsen & Richtler 
2000 and Larsen 1999). Very similar entities have also been 
found in interacting galaxies. Perhaps the best example of 
these is the Antennae with a collection of young star clus- 
ters with a median effective radius of 4 pc and ages of 
about 7 Myr (Whitmore et al. 1999). These galaxies also 
show other much larger entities (with an outer radius of 
450 pc) which arc not addressed in this paper. This star- 
forming activity in which masses similar to the total gas 
content found in galactic giant molecular clouds are turned 
into stars, all in a very small volume (~ few pc) implies 
the rapid accumulation of matter before star formation 
and negative feedback affect the collapsing cloud. Here 
we suggest that the feedback from continuously created 
massive stars (M* > 10 M©) can sustain a fragmenting 
small radius standing shell, giving origin to concentrated 
massive stellar clusters, with a universal IMF. 

Formation of a bound stellar cluster from the super- 
shell expanding out of the gaseous cloud has been dis- 
cussed by Brown et al. (1991, 1995). In their concept 
stars are formed when the supershell has swept out the 
entire cloud and expanded beyond its original boundary. 
However, the physical mechanism responsible for the su- 
pershell fragmentation and formation of the second gen- 
eration of stars remains uncertain. Here we give a thor- 
ough discussion of the supershell gravitational stability 
and show that the standing shell that forms due to the 
balance between the ejected mass ram pressure and the 
central cluster gravitational pull may become gravitation- 
ally unstable allowing the continuous formation of the new 
generation of stars. 

Section 2 describes the stellar factory, its self- 
regulation and the physics that lead to its closure, once a 
massive compact cluster has formed. Section 3 deals with 
the spectrum of fragments (sizes and numbers) expected 
from the factory model and compares this result with ob- 
served values of the IMF. Self-contamination is analyzed 
in Sect. 4, where strong predictions on the metallicity of 
the resultant clusters are given. Section 5 deals with other 
observables such as the expected IR luminosity and the 
(Ha) equivalent width of the resultant clusters. Finally, 
our conclusions are drawn in Sect. 6. 



2. The star-forming factory 

Our model invokes pressure-bounded, self-gravitating, 
isothermal clouds, which may become gravitationally un- 
stable if sufficiently compressed (Ebert 1955; Bonner 
1956). The gravitational instability allows a large cloud 
(Mc ~ lO'^ - lO'^ Mq) to enter its isothermal (Tc ~ 100 
K) collapse phase (Larson 1969; Bodenheimer & Steigart 
1968; Foster & ChevaHer 1993; Elmegreen et al. 2000), 
thereby developing a density and velocity structure with 
the following characteristics: 1) A central region of con- 
stant density (the plateau) where the velocity increases 
linearly from zero km at the center, to a maximum 
value {vma,^ = 3.3 Cc where Cc ~1 km s~^ is the sound 
speed of the collapsing cloud) at the boundary. 2) A re- 
gion of increasing size and constant maximum velocity 
{vmax), where the density falls off as (the skirt). As 
the collapse proceeds and the density in region 1 becomes 
larger, the knee region in the density distribution, that 
separates zones (1) and (2), moves closer to the center of 
the configuration with an increasing speed. As the density 
in the plateau region (pp) increases, unstable fragments 
begin to form. These will first (as pp grows larger than 
10^2" g cm^'') have a Jeans mass similar to those of mas- 
sive stars: Mjeans(g) = 8.5 x 10^^ " p-O-^. 

Although the above equation does not take into ac- 
count turbulence and magnetic fields as in Mac Low & 
Klessen (2003), one can, as a first approach, assume that 
a first stellar generation with a sufficient number of mas- 
sive stars (M» =100 - 10 Mq) will form spontaneously in 
the central plateau region where the 3D converging flow 
may trigger their gravitational instability. The stellar frag- 
ments will detach from the flow and, on time-scales of 
the order of a few times 10'"' yr, will enter the main se- 
quence. From then onwards, through their winds and ter- 
minal supernova (SN) explosions, they will begin to have 
an important impact on the collapsing cloud. For this to 
happen however, massive stars ought to form in sufficient 
numbers as to jointly stop the infall at least in the most 
central regions of the plateau. Otherwise, individual stars, 
despite their mechanical energy input rate, will Tinavoid- 
ably be buried by the infalling cloud, delaying the impact 
of feedback until more massive stars form. We thus as- 
sume that the first generation of massive stars is able to 
regiilate itself by displacing and storing the high-density 
matter left over from star formation into a cool expand- 
ing shell, thereby limiting the number of sources in the 
first stellar generation. One can show that, given the high 
densities (n > 10^ cm~^) attained both in the wind and 
in the plateau region and the size of the latter (~ 1 pc) a 
mass ~ 10^ Mq is available for a first generation of stars 
(with masses > 10 Mq). The shell will be driven by the 
momentum injected by the central wind sources, and the 
ionization front will rapidly become trapped within the 
expanding layer (see Sect. 2.1). Thus the two most dis- 
ruptive agents thought to interrupt (Hoyle 1953) or even 
stop (Cox 1983; Larson 1987; Franco et al. 1997) the for- 
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mation of stars are kept well under control by the infalling 
cloud. The large densities also promote a rapid radiative 
cooling within the shell, allowing low temperatures (~ 10 
K). The expanding layer is at all times confronted with 
the increasing density and larger velocity of the matter 
in the unperturbed plateau region, and it soon becomes 
ram-pressure confined. That is, it will soon happen that 
the central wind ram-pressure (pw^^w) exactly balances the 
infalling cloud plateau ram-pressure. From then onwards 
the shell of swept-up matter is forced to recede towards 
the stars, given the increasing density and velocity of the 
undisturbed collapsing plateau region, causing eventually 
the burial of the first stars. A more interesting situation 
arises if the number of sources in the first stellar genera- 
tion is such that the final position of the layer of shocked 
matter is close to the knee of the density distribution 
(Rk), where both the infalling gas density (pk) and veloc- 
ity (I'max) attain their maximum values. There, near R]^, 
the wind ram-pressure (pwV^) will exactly balance the in- 
falling cloud ram-pressure (pki'max) when the mechanical 
luminosity (L*) of the first spontaneous stellar generation 
approaches the critical value 



Lcrit = 2ttRIp]^v^vI^. 



(1) 



Star formation will suddenly stop as all plateau matter 
left over from star formation is now locked in the standing 
shell. Once the shell acquires this standing location, it will 
be able to process the infalling cloud mass. In this way, 
the mechanical energy deposited by the first generation 
of massive stars favors the accumulation of infalling cloud 
mass in the standing shell, which becomes gravitationally 
unstable. 

In our steady-state model everything happens at the 
same time. Gravitationally bound fragments continuously 
form in the unstable shell (at R = i?k) and then, due to 
their negligible cross-section, freely fall towards the center 
of the configuration as they evolve into stars. The larger 
number of sources continuously enhances the mechanical 
luminosity of the forming cluster and with it the amount of 
mass returned as a wind into the shell (Mw). At the same 
time, the continuous fragmentation of the shell and the 
infall of the resultant fragments acts as a source of mass 
in the most central region of the collapsing cloud, and 
this rapidly modifies the balance previously established 
between the wind and the infalling gas ram pressures. 
Indeed the ram pressure exerted by the wind sources, in 
order to keep the shell at its standing location, will now 
have to balance not only the infalling gas ram pressure 
but also the gravitational force exerted on the shell by the 
increasing mass of the central star cluster: 



AttRIp^vI 



AttRIp^v"^. 



(2) 



In the steady-state regime considered here, the cen- 
tral star cluster mass, Mgc, and the shell mass, A/gh, are 
Msc = 47ri?^pkWmaxt and M^h = 47ri?^Esh, where Esh is 
the shell surface density and t is the evolutionary time. 



Comparing the first and the second terms on the right- 
hand side of Eq. ((SJ one can show that very soon, after 
t > Vmax/ (47rGSsh) ~ 10"* — lO^yr, the infall ram pressure 
becomes negligible compared to the gravitational pull pro- 
vided by the forming cluster. Thus the shell becomes grav- 
itationally bound with its mechanical equilibrium simply 
given by the equation 



2 2 GM,hM„ 
47ri?i,pwi'w = 



i?2 



(3) 



One can show that despite the increasing effective gravity, 
the shell remains stable against Rayleigh - Taylor modes, 
because the density of the shocked wind is larger than 
that of the shocked infalling cloud. Also, nonlinear thin 
shell instabilities (see Vishniac 1994) and their induced 
pressure perturbations will be overcome by gravitational 
forces. 

From the new equilibrium condition (Eq. |2l , one can 
find how the ejected mass density grows with time at the 
standing radius r = i?k 



-t : 



max 



(4) 



and thus determine the mechanical energy input rate, Loq 
= ^Mw(i)'^wj required to support the shell in its equilib- 
rium state: 



in 



1 



STT^GSshPk-RfcWwWmaxi. 



(5) 



i.e., to support the shell against the gravitational pull ex- 
erted by the forming central star cluster, the mechanical 
luminosity would have to grow linearly with time. 

A second constraint on the mechanical luminosity 
arises from a consideration of the star formation rate 
(SFR). This is defined by the sum of the two sources of 
mass at the shell: the rate at which the collapsing cloud 
is processed by the shell (Min), which is a constant, and 
that rate at which mass is ejected by the star cluster (M^), 
which increases linearly with time. Thus, 



SFR{t) = Afin + Mw = 47ri?^pk^'max 1 



4^GI]sh 



(6) 



is a function that increases also linearly with time. Such 
star formation rate defines the energy deposition (£sc)j ex- 
pected from the growing central star cluster. Comparing 
Leq (Eq- HI) with isc, one can then find the value of the 
shell surface density Egh- -^sc is to be derived from star- 
burst synthesis models (Leitherer & Heckman 1995, Mas- 
Hesse & Kunth 1991, Sihch et al. 2002) taking into con- 
sideration the SFR prescribed by (Eq. |HJ and an assumed 
metallicity ISM of the host galaxy. 

Figure 1 shows that the energy input rate derived 
independently from the starburst synthesis models is in 
reasonable agreement with the equilibrium value (Eq. 
0) over a considerable span of time, particularly when 
Ss/i ~ 0.7 — 1.1 g cm^^. In such cases both mechanical 
energy input rates agree to within less than a factor of 
two - three over almost 25-30 Myr. 
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Fig. 1. Mechanical energy requirements. The figure shows 
the ratio Lsc (Eq. 6) over Lcq (Eq. 5). The Lsc values re- 
sult from starburst synthesis models that assume a SFR(t) 
as given by Eq. 6, for different values of Esh = 0.5 (dot- 
ted lines), 0.7 (solid lines) and 1.1 g cm"'^ (dashed lines). 
All models also assume upper and lower mass limits equal 
to 100 Mq and 1 Mq and a slope of -2.25 for the high 
mass end (as derived in Sect. 3). Lgc assumes mechanical 
energy input rate from a continuous star formation rate 
(Leitherer & Heckman 1995) assuming a galaxy metallic- 
ity value — 0.25 Zq. 



A star-forming factory then results from a profound 
self-regulation that accounts for the mass continuously 
added to the forming cluster, as well as for the mechani- 
cal energy that results from this further addition of mass 
and its transformation into stars. Self-regulation keeps the 
shell at its standing location and thus with the same frag- 
menting properties, while the forming cluster remains hid- 
den behind the shell and the collapsing cloud. 



2.1. The negative feedback caused by photoionization 

The increasing star formation rate also leads to a rapidly 
growing number of ionizing photons (Nsc)- This has also 
been calculated (see Fig. 2a) using the starburst synthe- 
sis model (Silich et al. 2002) under the assumption of 
a linearly increasing star formation rate, as prescribed 
by relation 10. Clearly, the ionizing radiation may up- 
set the shell fragmenting properties by simply changing, 
through photoionization, the temperature of the swept- 
up gas. Following Comeron (1997) and Tenorio-Tagle et 
al. (1999) we have derived the impact that such an in- 
creasing ionizing photon flux has on the collapsing shell. 

The shell density has been calculated from the momen- 
tum balance: 



where pw as given by Eq. is a function of time, and 
Csh is the sound speed of the shell. The full thickness of 
the shell is 



(8) 



where /in and iigh are the mean mass per particle and the 
shell number density, respectively. Thus the number of 
photons required for a complete ionization of the shell is: 



47r 



(9) 



Psh = Pw{Rk) — 
Csh 



(7) 



where ap = 2.59 x 10^^"^ cm"^ s^^ is the recombination 
coefficient to all levels but the ground state (Osterbrock, 
1989). A comparison of the critical number of photons re- 
quired for a complete ionization of the shell, with that of 
ionizing photons (A^sc) emitted by the growing central star 
cluster (Nsc/Ncrit) as function of evolutionary time (i) is 
shown in Fig. 2b. Throughout the evolution, this always 
remains ^ 1. This implies is that the ionization front, de- 
spite the increasing number of photons, is trapped within 
the shell structure and furthermore, is only able to photo- 
ionize a narrow inner section, the inner skin, of the gravita- 
tionally unstable shell. It is worth noticing that a fraction 
of the UV flux will be absorbed within a free wind region, 
reducing further the number of the UV photons reaching 
the shell per unit area and per unit time. Therefore the 
two possible negative feedback mechanisms, the ionizing 
radiation and the mechanical energy of the forming clus- 
ter, remain at all times under control by the star- forming 
factory. 

2.2. The size of the resultant clusters 

The factory stops operating either because small clouds 
(Mc < 10^ Mq) are rapidly processed by the standing shell 
or, in the case of larger clouds, because the mechanical 
energy input rate implied by the SFR condition (Lgc), 
here derived using starburst synthesis models, begins to 
largely exceed the luminosity required (icq) to keep the 
shell in its standing location. This latter possibility arises 
after 25 Myr of evolution (see Fig. 1), when the luminosity 
generated by the increasing SFR begins to dominate the 
equilibrium condition, although not by more than a factor 
of 2. After this time the shell will lose its standing location 
and be disrupted as it accelerates into the skirt of the 
remaining cloud. Thus, the size of the largest resultant 
clusters in our factory model, is restricted to a few IQ^Mq, 
the amount of cloud mass that can be processed by the 
standing shell within this time interval. 

Consequently, the half-light radius of the resultant 
cluster corresponds to the radius of the standing shell, 
and to the fraction of mass of the remaining cloud which 
is accelerated and removed from the original cloud when 
the mechanical energy input rate Lsc surpasses the equilib- 
rium condition. Assuming that 1/2 of the original cluster 
mass has been removed, the radius of the cluster doubles, 
reaching the value of a few pc. The final size of the cluster 
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Fig. 2. The effects of photoionization. a) The ionizing ra- 
diation produced by the forming cluster (Nsc) as a func- 
tion of time as derived from the starburst synthesis models 
under the same assumptions as those used in Fig. 1 and 
a Esh = 0.7 g cm~^; b) a comparison of the ratio of the 
photon flux (Nsc) and the number required to fully ionized 
the star forming shell (A^crit), as a function of time. 

depends on its subsequent internal evolution and on the 
environment in the home galaxy. 

Let us assume a cloiid of 10*^ Mq that has devel- 
oped the plateau-knce-skirt structure during its isother- 
mal {Tc = 100 K) collapse phase. By the time the knee 
reaches a radius of 2 pc, pp would be ^ 10^^" g cm^"^, and 
massive stars (M* < 100 Mq) will begin to appear at the 
center of the collapsing configuration. These will store the 
surrounding gas into an expanding shell, limiting the num- 
ber of sources in the first stellar generation. The expanding 
shell becomes ram-pressure confined within the plateau. 
However, given the increasing density in the unperturbed 
plateau region, the shell must recede until it finds the knee 
of the density distribution, where both density and veloc- 
ity of the incoming gas remain at a constant value. Here 
we assume a final shell standing at a distance of 1 pc. The 
density pk would then be « 6.68 x 10~^° g cm~^, and the 
infalling velocity equals 3.9 km s~^. At this time the outer 
radius of the 10*^ M© cloud is R^ax ~ Mc/(47rpki?k) ~ 90 
pc. Assuming = 10^ cm s~^, the mechanical luminosity 
derived from Eq. (1) is i w 5.8 x 10^^ erg s~^. 
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Fig. 3. The dispersion relation Lu{ri) for Csh correspond- 
ing to Tsh = 10 K, i?k = 1 pc and E^h = 0.7 g cm^^. 
Vmax = 503, corresponds to the wavelength A = 0.01 pc, 
gravitationally unstable mass m = TrSghA^ = I.6M0, and 
the e-folding time of the growth of perturbations of 5 x 10"^ 
yr. The largest and smallest masses form at 77 = 62 and 
1005 corresponding to the mass range 0.39 - 150 Mq. The 
upper axis indicates the (non-linear) mass scale, and the 
arrow the position of a 100 M0 fragment. 

In the steady-state regime, the infalling mass ends up 
being transformed into stars at the rate at which it is 
processed by the standing shell Mjn = 47ri?^pkWmax ~ 
4.4 x 10~^ Mq yr~^ An originally 10^ Mq cloud would 
then be processed in about 24 Myr, while the mechanical 
luminosity of the forming cluster increases almost linearly 
with time up to 6.3x10'*" erg s~^. 

Given the self-similar solution of the isothermal col- 
lapse phase, all collapsing clouds are processed in a very 
similar manner, all with a similar IMF, and all of them 
causing a similar SFR rapidly increasing with time. Thus 
the only limitation on the mass of the resultant clusters 
is set by the time during which they are able to enhance 
their mechanical luminosity with time, to keep the frag- 
menting shell at its standing location. Note also that for a 
given plateau density, collapsing clouds present the same 
size plateau, regardless of the mass of the collapsing cloud 
(Mc). Thus, as the radius of the star-forming shell is in- 
dependent of the collapsing cloud mass, the mass-radius 
relation observed for GMCs is not reflected in the stellar 
clusters, as pointed out by Ashman & Zepf (2001). 

3. The mass spectrum of fragments 

Mass accumulation leads to the gravitational instabil- 
ity of the standing shell with a well-defined mass and 
number of resultant fragments. The dispersion relation 

for gravitational instability of an expanding shell of ra- 
dius R is (Elmegreen 1994; Wiinsch & Palous 2001): 

- = -f + (f - ^ + '-^f" , where S., is 
the unperturbed surface density of the shell, Csh its sound 
speed and G is the gravitational constant. The condition 
for instability demands w to be real and positive. The 
wavenumber r] is related to the perturbation wavelength 
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X hy T] = 2ttRX^^, and the e-folding time of the pertur- 
bation growth is uj^^. For a standing shell configuration 
(with _R = 0) this reduces to 



i?2 



27rG'i:sh?? 
R 



(10) 



The dispersion relation 0^(77) is shown in Fig. 3. The 
fastest growing mode occurs for Tymax — 



corre- 



sponding to LJ-n 



ttGS. 



^. Fragments of mass 



(11) 



form with a frequency proportional to the growth rate, 
which is given by lo, and is proportional to the shell sur- 
face available to accommodate fragments of a given wave- 
length, -p- = Consequently the number of fragments 
AiV formed per unit time, out of the standing shell, with 
the wavenumber from the interval (77, 77 4- A77) is 

/ 277 77^ <:\uj\ ^ 

The A77 interval, can be expressed in terms of the corre- 
sponding A777 by means of Eq. ((m) . And thus, Eq. ((TT)l 
and p2ll imply the mass spectrum of fragments, or define 
the initial mass function (m) = AN/ Am. 



'2777-9/4 



1/2 



(-7rl/2c2j^,„-l/2+GI]l/2)l/2' 



(13) 



where Q is the normalization factor fixed by the star for- 
mation rate (SFR). Note that a negative Ar77, should be 
used for positive At;. Thus, mass accumulation leads to the 
gravitational instability of the standing shell with a well 
defined mass and number of resultant fragments. From 
the dispersion relation of the linearized analysis of the 
hydrodynamical equations on the surface of the standing 
shell, the mass spectrum of gravitationally bound frag- 
ments presents a slope equal to -2.25 for massive objects. 
The distribution flattens in the neighborhood of rjmax and 

ire* 

peaks at m — ^Q'i"^ ^ ■ The minimum mass, obtained from 
the condition that to is positive, lies at 2ri^i^^. These results 
are in good agreement with the stellar mass distribution 
(see Fig. 4) inferred for star clusters (Moffat 1997; Hunter 
et al. 1997; Wyse 1997; Piotto & ZoccaU 1999; Paresce & 
De Marchi 2000) and for the solar neighborhood (Salpeter 
1955; Scalo 1986; Binney & Merrifield 1998; Kroupa 2001, 
2002). 

4. Self-contamination 

After ~ 3 Myr of evolution, the mass ejected from the star 
cluster (Mw) becomes rapidly contaminated by the super- 
nova explosions. We assume that on reaching the shell the 




m(Me,) 

Fig. 4. The IMF - ^ (m) - as given by Binney & Merrifield 
(1998), is compared with our results from Eq. (13) normal- 
ized to the total cluster mass (solid line) . The comparison 
assumes a 10^ Mq cloud, used below as an example, fully 
processed into stars during a time span of 24 Myr (see 
text) 

metal-enriched matter is immediately mixed with the gas 
coming from the collapsing cloud (Mm) which, depending 
on the host galaxy, may present a low metallicity. This 
leads to a continuously increasing abundance for every 
new generation of stars resulting from the fragmenting 
shell. 

We assume that the gas ejected by stellar winds and 
SNe includes all the newly synthesized metals as well as 
metals distributed throughout the stellar hydrogen en- 
velopes of the progenitor. The total mass ejection rate 
and the rate of metal ejection Mmet for a central star 
cluster with a power-law initial mass function then are 



M, M{t + At) 



2-Q 



M{t) 



2-Q 



^ At 



M2-" _ A//"2"" 



(14) 



(«-2) 



E 



Ml 

At 



M{t) 



M(t+At} 



(15) 



where Mi is the mass of every new generation of stars, t is 
the evolutionary time. At is the time step, Afiow and M^p 
are the lower and upper cut-off masses, respectively, M(t) 
is the mass of the stars exploding after an evolutionary 
time t (see Silich et al. 2002). 

The evolution of the iron abundance of the forming 
star cluster (which can be compared with the available 
metallicities of galactic super-star clusters) can be ap- 
proximated by means of an iron yield Ype- Here we use 
the Thielemann et al. (1992) model and extrapolate their 
results as constant yields for low (< 13 Mq) and high 
(> 25 Mq) mass stars (see for details Silich et al. 2001). 

We then adopt a helium normal abundance (one he- 
lium atom for every ten hydrogen atoms) and calculate 
the number of hydrogen and iron atoms mixing inside a 
shell at every time step 



iVu = 



1.1 fir 



(16) 
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ApeniH ' 



(17) 



where /Zn = is the mean mass per particle, ^Fe = 

55.4 is the atomic number of iron in hydrogen mass units, 
and mn is the hydrogen mass. The total mass input rate 
Mtot) as well as the iron mass input rate Mpe, includes 
both the matter ejected as a function of time by the star 
cluster and the mass coming from the collapsing cloud: 

Mtot(t) = Mw+Min, (18) 
MFe(i) = AfFe.w + MFe,i„. (19) 

The iron abundance of every new generation of stars is 



Fe 

IT 



log 



Fe 



log 



Fe 



(20) 



where the Solar iron abundance is log (^j^^ +12 = 7.448 

(Holweger, 2001). Fig. 5 presents the calculated iron abun- 
dance of the fragmenting shell, as a function of time, for 
galaxies with different ISM metallicities. As the clusters 
form from the ISM gas as that is contaminated by the 

supernova products from former stellar generations, the 
resultant metal spread in a cluster depends strongly on its 
final mass (upper axis in Fig. 5). In this way, clusters with 
a final mass smaller than 2.7 x 10^ M©, processed within 
3 Myr, will display the metal abundance that would reflect 
that of the ISM at the moment of formation. More massive 
clusters however, as their formation time may exceed the 
supernova era 3 Myr) from former stellar generations, 
will display a large metal spread in their sources. A cluster 
of 10^ M0 (indicated by an arrow in Fig. 5), requires 24 
Myr to complete its formation. If such a cluster forms in 
a low metallicity (Zism = O.OIZ©) galaxy for example, it 
will present stars with different [Fc/H] abundance within 
the range —2 < [Fe/H] < —0.13, while equally massive 
clusters forming out of a more metal-rich ISM will show 
a correspondingly smaller relative metallicity spread (see 
Figure 5). Note however that the spread caused during for- 
mation of the clusters depends strongly on the assumed 
metal yields. 

5. Further observational properties 

Since the star-forming shell is dense, let us assume that 
there is enough dust to obscure the optical emission of the 
forming star cluster and transform it into the IR. The ex- 
pected IR nebular luminosity would then result from the 
stellar ionizing flux and from the degraded cluster me- 
chanical luminosity, thermalized at the reverse shock and 
then radiated away via effective gas cooling. Therefore the 
luminosity of the IR nebula associated for example to a 
~ 10^ Mq cluster may reach ~ lO^L©. 

After formation (see Sect. 2.2) the newly formed star 
cluster will become visible in the optical line emission re- 
gion. At such a time, the luminosities of a 10^ Mq and 
10^ Mq star clusters would be Lna ~ 1-4 x 10^° erg s~^ 




Time (Myr) 

Fig. 5. The [Fe/H] metallicity of the resultant clusters. 
The star forming shell metallicity as function of time for 
galaxies with different initial metal abundance: .^ism = 
0.01^0 (dashed hne) Zism = O.l^o (Sohd line) and 
■^iSM = Zq (dotted line). The resultant clusters, depend- 
ing on their mass (upper axis), will show a metal spread 
that would strongly depend on the ISM original metallic- 
ity. 

and Lua ~ 9.4 x 10^° erg s~^, respectively, if all ioniz- 
ing photons are trapped within the surrounding gas. The 
Lyman continuum rate of the 10^ M© cluster is around 
6.9 x 10^^ s~^, comparable to the Lyman continuum rate 
(3x 10^^ s~^) derived from the radio emission flux of super- 
nebula NGC 5253 (Gorjian et al. 2001). 

The continuous creation of stars in the star-forming 
factory model also leaves a foot-print on its Hq, equiva- 
lent width. Figure 6a displays the predicted Ha equivalent 
width as a function of the star cluster mass, at the mo- 
ment at which formation is completed and the cluster be- 
comes visible for the first time (Wop^j]). The continuous 
creation of stars during the formation phase leads to low- 
mass clusters with initial equivalent widths (M o[Ha]) 
larger than those arising from more massive clusters (see 
Fig. 6a). 

After formation the number of UV photons drops 
rapidly and the equivalent width accordingly decreases 
(Fig. 6b). The evolutionary plot implies an HII region life- 
time which is slightly shorter than 10 Myr after cluster 
formation, with an initially smaller equivalent width 
for larger masses of the resultant star clusters. 

6. Conclusions 

The formation of compact and massive stellar clusters is 
naturally explained in the framework of the star-forming 
factory. The shell that stores the collapsing cloud is able 
to maintain its standing location and fragmenting prop- 
erties, thanks to the larger energy input rate that results 
from the increasing number of massive stars in the central 
region. With this energy input rate, the forming cluster is 
able to balance the increasing gravitational pull exerted 
on the standing shell by the continuously growing mass 
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Fig. 6. The He, equivalent width, a) The Hq equivalent 

width Wq predicted at the moment when formation is 
completed and the clusters becomes visible, as function 
of the exciting star cluster mass. The synthesis models as- 
sumed the same upper and lower stellar mass limits, as 
well as the increasing star formation rate and the slope 
(a = -2.25) as in Fig. 1. b) The post-formation time evo- 
lution of the Ha equivalent width for star clusters with a 
mass of 10^ Mr and 10*^ M«. 



host galaxy ISM metallicity at the moment of formation, 
to values similar to solar metallicity (see Fig. 5). 

Upon formation, the more massive the resultant clus- 
ter, the smaller its initial H„ equivalent width. This may 
differ by factors of three between clusters with a final stel- 
lar mass of 10"'' and 10® M0. Large IR luminosities (> 10® 
Lq are predicted for factories leading to large (> 10® M©) 
super-star clusters. 

A definite prediction of the factory model, applicable 
in particular to high mass clusters (say ~ 10® Mq), is 
the possibility of finding, after formation, a mixture of 
stellar populations with different ages and metallicities. 
For example, within a massive cluster the model predicts 
the co-existence of O stars and supergiants evolving at the 
same time as WR stars, while other stars may explode as 
SN. All of this is the result of the evolution of consecutive 
generations of stars born at different stages during the 
formation of the cluster. 
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in stars within the central region. Self-regulation in the 
factory can only be sustained for up to 20 - 25 Myr, while 
the mechanical energy input rate grows to satisfy the re- 
quired energy to keep the shell in its standing location. 
Afterwards, the energy input rate grows and exceeds the 
equilibrium condition, causing the acceleration of the shell 
and with it its disruption as it moves into the skirt of the 
collapsing cloud. After that time, a few times 10® Mq 
would have been converted into stars, all following a sim- 
ilar IMF. Smaller proto-cluster clouds are fully processed 
on shorter time-scales. 

Small clusters (< 2.7 xlO^ Mq) would present the 
metallicity identical to the metallicity of the ISM out of 
which they formed. More massive clusters however, will 
carry the footprint of self contamination, produced by the 
supernovae products from former stellar generations used 
to support the factory at work. Massive clusters will then 
present a metallicity spread that would range from their 
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